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Abstract We present different methods for symbolic computer algebra compu- 
tations in higher dimensional (> 9) Clifford algebras using the CLIFFORD and 
Bigebra packages for Maple®. This is achieved using graded tensor decomposi- 
tions, periodicity theorems and matrix spinor representations over Clifford numbers. 
We show how to code the graded algebra isomorphisms and the main involutions, 
and we provide some benchmarks. 



1 Introduction 

Clifford algebras are used in several areas of mathematics, physics, and engineering. 
Since computing power has increased tremendously, in terms of available memory 
and actual processing speed, practical symbolic computations, say on a laptop fin- 
ishing within minutes, are by now possible for Clifford algebras over vectors spaces 
of dimensions higher than 8. When CLIFFORD was designed, way back in the early 
90's m, such computations where impossible. However, the design of CLIFFORD, 
restricting the vector space dimension to less than or equal to 9, incorporated from 
the beginning the idea that using mod-8 and other periodicity isomorphisms of real 
Clifford algebras allows nevertheless computations in higher dimensional algebras 
as well. Using periodicity theorems also 'groups' Clifford algebras as in a periodic 
table, the spinorial chess board ifTOl . Hence, implementing the periodicity in sym- 
bolic computations -as described in this paper- makes use of these intrinsic algebra 
features. Recent applications in engineering, for example when modeling geomet- 
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ric transformations in robotics, rely on real Clifford algebras like C%2 Q and thus 
enforce the need for using the periodicity in computations with CLIFFORD. 

In this note, we will describe three ways showing how CLIFFORD deals with 
(real) Clifford algebras of higher dimensions (> 9). One way is to use Bigebra, 
an extension of CLIFFORD, to utilize graded tensor products of Clifford algebras. 
Another method utilizes ungraded tensor products and periodicity isomorphisms 
without a need for introducing spinorial bases. The third method uses, on top of the 
periodicity theorems, a spinor representation of one factor of the tensor decompo- 
sition, hence computing in Clifford algebra valued matrix rings. For example, from 
Cip+\,q+\ — Clp,q®Ci\^\ ~ Mat(2, C£ p ,q) we see that computing say in Cig^ can 
be done with the 2x2 matrices over Cijj. CLIFFORD supports computing with 
matrices having Clifford entries. Similar computations still make sense, e.g., for 
conformal symmetries using Vahlen matrices Mat(2,C£3j) ~ Q4.2. 

After recalling our basic notations and the periodicity theorems, we start explain- 
ing how we use Bigebra to compute with graded and ungraded tensor products 
of Clifford algebras. Then we proceed to the third matrix-based method indicated 
above. As this method relies on spinor representations of real Clifford algebras, one 
encounters not only real, but also complex, quaternionic, double real, and double 
quaternionic spinor representations. For the sake of simplicity and space, we will 
just deal with real representations of simple Clifford algebras. The periodicity theo- 
rems single out the signature cases where a graded algebra isomorphism is available 
onto an ungraded tensor decomposition and hence, introducing spinor bases, matrix 
tensor products can be employed. 



2 Tensor product decompositions and periodicity for C^-algebras 

We study Clifford algebras over finite (real) vector spaces. While these algebras can 
be defined by a universal property, for actual computations in a CAS (Computer 
Algebra System) we use generators and relations. 



2.1 Basic notations, quadratic and bilinear forms 

Let V be a finite dimensional real (or complex) vector space with scalar multipli- 
cation R x V — > V :: (A,v) A v. A quadratic form is a map Q : V — > R such that 
Q(Xv) = X 2 Q(v), with associated polar bilinear form B : V x V — > R derived from 
Q and defined as 2B(v, w) = Q(v + w) - Q(v) - Q(w) (hence Q(v) = B(v, v)). Q is 
called non-degenerate if Q{v) = implies v = (v G V,Vw G V,B(w,v) = implies 
v = 0). An isomorphism V ~ R" = ©"R defines a set of generators for V from the in- 
jections of the direct sum : R — > R" :: 1 h-> e^ ~ v^, that is a basis for V. Sylvester's 
theorem states that there exists a basis for V such that the quadratic form Q is di- 
agonal with entries ±1,0 in the real case (0 only when Q is degenerate; just +l's 
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in the non-degenerate complex case). Under these isomorphisms, the real quadratic 
space (Q,V) is isomorphic to a space R /, - t/ = (R n ,5 M ) with a diagonal quadratic 
(polarized) form 8 p . q with p ones and q minus ones. (p,q) is called the signature 
of Q and p + q = n = dim V. Furthermore we have generators e k of W' q such that 
Q(e k ) = +1 for 1 < k < p and Q(e k ) = -1 for p < k < p + q. 



2.2 Grassmann algebra, Z- and Z2-gradings, main involutions 

We can associate functorially to any (say finite, real) vector space V the Grassmann 
algebra of antisymmetric tensors, /\ : V — > A^ = ®"=o A'^- Using the linearly or- 
dered basis {e k } of V (e k < e/ if k < I), we get a basis 

{e,} = {e,-j A e,- 2 A • • • A e,- ; |7 = (/'i ,i 2 , . . . ,//), ;'i < ; 2 < ••• < //, ij G {l,...,n}} 

for the space f\V of antisymmetric tensors. We can extend the order on V to the 
inverse lexicographic order on /\V. We associate to e/ a degree \I\ (number of gen- 
erators e k in ej) and a grade (or parity) as |/| mod 2. The Grassmann algebra is 
Z-graded w.r.t. the degree, that is f\'V f\ J V C f\' +J V. This can be restricted to a 
Z 2 -grading /\V = Ao^® Ai^ w.r.t. the grade, hence decomposing into even and 
odd parts (Ao^ is a subalgebra, Ai V is a Ao^~ mo dule). The subspaces A' V have 
(") many basis elements and hence A V has dimension 2" = (") . Homogeneous 
elements vi A ■■■ A v, of A' ^ are called extensors (or blades). A general element 
X 6 A^ is an a gg re g ate ^ = L x / e / w i tn rea l coefficients xj. A presentation using 
generators and relations is given as 

f\V := (e k | e k ei = -e,e^, i,k G {l,...,n}). 

The Grassmann algebra comes with two main involutions. The grade involution 
extends the map — : V — > V :: v t— > — v (additive inverse in V) to (") on f\V. On 
generators this reads as ej = (— 1) ' 7 'e/ and it extends by linearity. For example, e\ = 
2.3 := Ae? Ae3 = (— l) 3 <?i. 2 .3 = — e l,2,3 while e\ >2 '■= <?i Ae 2 = e i.2- 

A second involution makes use of the opposite algebra. Let (A,m.A) with : 
AxA—tA :: (a,b) M> ab be an algebra. The opposite algebra A" p = (A,m'^ p ) is 
the same vector space A with the multiplication m op : A x A — > A :: (<z,Z>) i-> £>a. 
The opposite wedge product A" p is given as m A op v = v Am (no signs), producing 
the opposite Grassmann algebra /\" p V. The reversion involution on A V extends the 
identity map 1 : V —> V to A V in such a way that rev o A" p = Ao (rev ® rev). In other 
words, rev : /\" p V — > A V is a Grassmann algebra isomorphism. As rev is invertible, 
this can be used to define A op . On basis elements we have e/ = (— l^KM -1 )/ 2 ^ 
with the usual notation rev(u) = m. 

Let V* = [V,R] be the dual vector space (linear forms). One finds that A^* — 
(A V)* is the dual Grassmann algebra, and we can define a pairing 
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'det{e*{e k )) : when |/| = |/|; 
0, when |/| ^ |/|. 



The interior product J is defined as the adjoint w.r.t. the duality pairing of multipli- 
cation in f\V*, tnat is (u | vJw) := (v* Au \ w). Let x,y 6 V and u,v,w <E /\V then 
J is defined operationally by the rules (Chevalley construction) 

xJy=(x\y), xJ (uAv) = (xJu) Av + «A(xJv), (u A v) J w = u J (v J w) . (1) 

Using a nondegenerate bilinear form B as duality, one obtains x J y = B(x,y). 



2.3 (Real) Clifford algebras 

Given a quadratic space (V, Q) we can functorially associate to it the universal Clif- 
ford algebra C£(V,Q) 03). Given the isomorphism (V,Q) ~ R p ' ? this provides the 
Clifford algebra C£(M.",8 P ^) (« = p+q) also often denoted as R p<? . As vector 
spaces one has W'-' 1 = /\W, hence we have a Grassmann basis {e/} spanning the 
vector space underlying the Clifford algebra. The Clifford product can be imple- 
mented in various ways. A Clifford algebra presentation reads 

C^P,q '■= i e i I ejej + ejej = 0,i^= j; ef = 1 if 1 < i < p, otherwise ef = — 1} 

where the Clifford multiplication is juxtaposition and i,j G {1, . . . ,«}. Another way 
to define the Clifford multiplication uses the interior multiplication and the Clifford 
map y. Let x S V, u S /\ V then y x u — xu =xJu + xhu and extend by Q} with 
x Jy = B(x,y) and linearity. Defining CI as a quotient of the tensor algebra it can 
be seen that it is not graded, but only ^-graded, CI = CIq ®C£i with C£q the even 
sub-Clifford algebra and Cl\ a C^o-module. In terms of the Grassmann Z-degrees, 
it is a filtered algebra C&CV C ^C^. The quadratic space (V, Q) decomposes 

as (V,Q) = {Vi+V2,Qi -LQi) with restricted quadratic forms Q\ =Q\y- l on V\ and 
Qi = Q\v 2 on Yi- The Grassmann algebra functor is exponential, that is, 

/\(V l +V 2 )=/\(V l )®/\(V 2 ) 

(the graded tensor product (g) is defined below and we use equality for categorical 
isomorphisms). Similarly we get for Clifford algebras 

C£(V l +V 2 ,Q l ±Q 2 )=a(V 1 ,Q l )&a(V 2 ,Q 2 ), 

and it is this decomposition which will be used below to compute in CLIFFORD 
in dimensions > 9. For Clifford algebras with non-symmetric bilinear forms such a 
decomposition is in general not direct, see lfl3l . 
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2.4 Tensor products of (graded) algebras 

Let (A , niA ) and (Z?,tob) be K-algebras. That is, we have right and left scalar multi- 
plications p : A x K — > A and X : IK x B — > B. This allows to define the tensor product 
A %B as a universal object via a co-equalizer c of two arrows pol and 1 o X 

AxKxB^-^rAxB — C -+A® K B , 

which produces the common relations such as aX <g) b = a <E) Xb, multilinearity etc. 
We have (vector space) injections i& : A — > A ® B and ig : B — !> A ® B, and want to 
transport the algebra structure too. We define the product algebra on A <EiB from the 
algebra structures on A and B as follows: m^B := (ffu ®mB){\ <8> o <8> 1), where 
(7 : A (g)B — > B(£)A :: (a,/?) n> is the switch of tensor factors. On elements: 

{a®b)(a' ®b') = (aa' ®bb'). (2) 

Note that we had to switch a' and b, assuming they commute. As we work in a 
^-graded setting A = Aq +A\, B = Bq +B\, this switch has to be replaced with 
a graded switch a : A® B — > B(§)A :: (a,b) h-> ( — l)l a ll fe l (b,a), (on homogeneous 
elements and extended by linearity). The graded multiplication m A ^ jB is then given 
by m A ^ B := (niA® ms)(l <E) (7® 1). Using this setup, the injections above become 
algebra isomorphisms i& : a >-t (a® 1) and (g : i-4 (1 <8 b). 

In the Grassmann algebra case, splitting the space V = V\ + V2 with n basis vec- 
tors e,- into two sets with, respectively, p (1 < i < p) and g (p < i < n) vectors, we 
get the maps e\ M- e; ® 1 (r < p) and ey i-)- 1 (§) e ; - (p < j < «). In the CAS compu- 
tations below we will standardize the indices, that is, we will reindex j h-> j — p so 
that ; £ { 1 , . . . , p} and 7 — p G { 1 , . . . « — p } . The graded tensor product ensures that 
we still have the desired anti-commutation relations 

(e,(§) 1)(1 <§)«,•) = (ei&ej) and (1 ® e 7 -)(e ; ® 1) = (-l) 1 ''^® e 7 -) . (3) 

In this way we can form graded tensor products of Clifford algebras Ci p ^® Cl r>s 
too, and that is what we aim for. Tensor products for matrix algebras are usually 
called Kronecker products, and are taken without the grading. Let A ~ Mat(«,R) 
and B ~ Mat(ra,R) be matrix algebras. Recall that matrices need a choice of basis 
giving matrices [«,■_,-] and [bjj] . The definition of the matrix tensor algebra A <8> B ~ 
Mat(« • m, R) includes a choice of how to form a basis for A ® B, which consists of 
elements E nm = eij^e^i, that is, a reindexing function [(i,j),(k,l)] M> (n,m). One 
way is to define a tensor [a, -j] ® = [a,-; • [»£/]] by inserting the matrix [Z^/] as 
blocks into the matrix [a.jj] ; another obvious choice would exchange the role of [a] 
and [b]. Category theory shows that the definition of the tensor algebra is unique 
up to a unique isomorphism depending on the particular choices. However, actual 
computations in a CAS need to be consistent and explicit in the choice of these 
isomorphisms. 
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2.5 Periodicity of Clifford algebras 

We have seen in the previous section that we can tensor Clifford algebras of any 
signature provided we employ the graded tensor product. 

Theorem 1. Let Cl pq , C£ r j be two real Clifford algebras, then 

C£p+r,q+s — Cip.q <§> Cd r ,s (4) 

( which does not even need nondegeneracy, or even symmetry, of the involved forms). 

This result will be used in section[3]to describe a general method using Bigebra to 
do practical symbolic CA computations in higher dimensional real Clifford algebras 
of any signature. 

Using matrices over Clifford numbers, like Mat(2,C£ P]9 ), needs considering un- 
graded tensors, as the matrix algebra tensor products are ungraded. Doing so em- 
ploys (graded) algebra isomorphisms described on the generators of the factor Clif- 
ford algebras inside the ambient Clifford algebra. This leads to the well-known pe- 
riodicity relations which are summarized in the following^. 

Theorem 2. For real Clifford algebras we have the following periodicity theorems 
and isomorphisms: 

1) Ciq+hp-l - Cl M ifp > 1 (see ED), 

2) Ct p , q+1 ~ a q , p+1 (see ED, 

3) Ct p , q ~ Cl p -4, q +4 ifp>4 (see MM), 

4) a p+4 . q ~ Ce p . q (E)CU,o ^ C£ M <g>Mat(2,H) (see ED), 

5) C£ p , q+4 ~ Ci p . q ®Cl 4 fi ^ C£ Pi9 <g>Mat(2,!I) (see ED), 

6) a p+hq+ i~a p , q ®a 1A ~a Ptq ®Mat(2,R) ~Mat(2,CX M ) (see El), 

7) C£ p , q+& ~ a M ®Mat(16,M) ~ Mat(16,CX M ) with Mat(2,H) <g>Mat(2,H) ~ 
Mat(8,R) (see EDDM3JD), 

8) C£ p+S . q ~ a Pi g®Mat(16,R) ~ Mat(16,CX M ) with Mat(2,H) <g> Mat(2,H) ~ 
Mat(8,R) (see EDDMUD)- 

Note that here all tensor products are ungraded and these structure results charac- 
terize the cases (signatures) where a graded isomorphism from the graded case to 
the ungraded case exists. The isomorphisms with matrix rings need spinor represen- 
tations and will be discussed briefly in subsection l2.6l 

By way of example we show a typical isomorphism on generators for the iso- 
morphism 6). 

Example 1. Let {e,} be the set of orthonormal generators for Cl p . q with ej = 1 for 
i e {1,. . . ,p} and ej = -1 for i e {p + 1, . . . ,p + q}, and let fj = l = -f\ be 
the orthogonal generators for Cl\\. Then, e,- 0/1/2, 1 <S>/i, 1 (8/2 form a set of 
generators for Ct p+hq+l . Indeed, since fif 2 = -f 2 fi and (/1/2) 2 = 1, for i,j £ 
{!,..., p + q} and i ^ j, we find the following familiar relations in Ci Piq ®Cl\y. 



For additional references on the periodicity theorems see [8 9 12 1 5 1 
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(«< ® /1/2) 2 = («< ® 0/1/2) = «? ® (/1/2) 2 



( 



1®1, if 1 <i<p; 
1(8)1, otherwise, 



(1 <g/i) 2 = 1 ®/f = 1 ® 1 = -1 ®/f = -(1 ®/ 2 )2, 

(ei<8>/i/ 2 )(e;®/i/ 2 ) + (e > ®/i/2)(ei®/i/2) = fae/ + e/e;) (/i/ 2 ) 2 = 0, 

(«* ®A/y (1 + (1 "8/1/2) - *< ® (-/i 2 /2) +«* ® (/i 2 / 2 ) = 0, 

(e«®/i/2)(l 8/2) + (1 ®fi){ei (8/1/2) = * (/i/ 2 2 ) + 4 (-/i/l) = 0, 

(l®/l)(l®/ 2 ) + (l®/ 2 )(l®/l) = l®(/l/ 2 +/ 2 /l)=0. 

For further details see, e.g., II10II171 . 

Theorem 3 ( II 1011 . Thm. 5.8). With the notation as above, let V 2 have dimension 
2k and let CO be the volume element in C£(V 2 ,Q 2 ) with CO 2 = X 7^ 0. There exists a 
vector space isomorphism between the module C£(V\ ©V^iGl -L Ql) and the module 
C^(Y\, jQi) ® CliViiQi) given on generators as (x,y) 1— > x ft) + 1 y, and there 
is a graded algebra isomorphism 



The involutions extend as (x0y) ~ x0y andrev(x®y) ~ rev(x) ® rev(y) ;/ |x| = 
mod 2 even andmv(x®y) ~ rev(x) rev(y) otherwise. Then all periodicity isomor- 
phisms in theorem\2\are special cases of this one^ 

To exemplify this, let (x,y) be any pair of generators with x G V\ and y G V 2 which 
upon the embedding V\ ©V 2 Gf?(Vi ©V 2 ,2i -L Q 2 ) we write as the sum x + y. 
Then, 



(x + y) 2 =x 2 + (x 3 ; + yx)+y 2 =(e 1 (x)+e 2 (y))l = (e 1 ±e 2 )(x,y) (6) 



due to the orthogonality of x and y. On the other hand, in the (ungraded) tensor 
product algebra in the right-hand-side of (|5) we find, as expected, 



due the anti-commutativity yft) = —coy for every y G V 2 assured by the even dimen- 
sion of V%. Furthermore, this last computation shows that the assumption ft)" = X ^ 
and the appearance of the factor i modifying Q\ in Cl(V\ , jQi )> are both necessary. 




(5) 



(x®ft)+l0y) 2 

= (x0ft))(x0ft)) + (x0ft))(l0y) + (l®y)(x®ft)) + (l0y)(l®y) 
= x 2 0ft) 2 +x® ft)y+x0yft)+ 1 0y 2 

= 02i(x) + C&Cy))(i®i) 



(7) 



2 The requirement (B 2 = A ^ is equivalent to (22 being non-degenerate. See also [US] p. 218]. 
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Since we are using Grassmann bases in all Clifford algebras, it is interesting to 
calculate the image of a Grassmann basis monomial, lets say of degree 2, in the 
product algebra on the right of (|5}- Let X\,X2 be orthogonal generators in V\ and 
let yi,y2 be orthogonal generators in V 2 - Then, the wedge product in the Clifford 
algebra Ci(Vi © V2,Gi -L Qi) is computed as expected 

(xi +yi) A(x 2 +y2) =x\ Ax 2 +xi Ayi+yi Ax 2 +y\ Ay 2 - (8) 

On the other hand, since (xi +yi) A (x 2 + ^2) = (*l + y\)( x 2 + yi), its image in 
C£(V\, tQi) ®Cf(V 2 ,<2 2 ) under the isomorphism ((5} is the following rather compli- 
cated element: 

(xix 2 ) ® CO 2 + (xi A 1) ® 0)y 2 + (1 Ax 2 ) ® yi to + (1 A 1) © (yiy 2 ) = 

(x] Ax 2 ) © Al +xi © 0)y 2 +x 2 (giyiCO + 1 ® (yi Ay 2 ). (9) 

The isomorphism in (O is given by the procedures bas2Tbas (from left to right) 
and its inverse Tbas2bas (from right to left). In the worksheets []6] we show both 
procedures as well as we verify the assertions regarding the involutions. 



2.6 Spinor representations, Clifford valued matrix representations 

A Clifford algebra is an abstract algebra, but we may want to realize it as a con- 
crete matrix algebra. It is, however, well known that matrix representations may 
be very inefficient for CAS purposes. The simplest representation is the (left) reg- 
ular representation, sending a £ A n> X a = »m(«,— ) £ End (A), the left multipli- 
cation operator by a. This representation is usually highly reducible. The small- 
est faithful representations of a Clifford algebra are given by spinor representa- 
tions. Algebraically, a spinor representation is given by a minimal left ideal which 
can be generated by left multiplication from a primitive idempotent fi = ff with 
flfkJi idempotents such that /; =f k +fi and f k f t = f,f k = 0. The vector space 
S := C£p, q fi is a spinor space, and it carries a faithful irreducible representation of 
C£ p ,q for simple algebras However, when C£ p , q is not simple, and in several signa- 
tures (p,q) this space is not really a vector space, but a module over K = fiCl p . q fi 
with K isomorphic to R, C, H, 2R = R © R, or 2H = H © H depending on the 
signature (p,q). The spinor bi-module ce^R carries a left Ci Piq and right K action 
(scalar product). Looking up tables of spinor representations |[T6l[T8l yields that 
C£ P} q is simple and has a real representation if p — q = 0,2 mod 8, distinguished 
by the fact that a normalized volume element © squares to (O 2 = + 1 if p — q = 
0, 1 mod 4 and ft) 2 = — 1 if p — q = 2, 3 mod 4. Avoiding non real K's, we find the 
isomorphisms C£ .o ~ M ~ Mat(l,R), {Ce 2 ,o,C£i A } ~ Mat(2,R), {C£ 3 .i,C£ 2 .2} ^ 
Mat(4,R), {Ce 4 ,2,a 3 3,a , 6 } ~ Mat(8,R), and {C4,o,« 5 3^4,4,«i, 7 ,C^ 8 } ~ 

3 For semi-simple Clifford algebras, we realize the spinor representation in S ffi S where S = G? M /; 
(see, e.g., fl6l ). 
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Mat(16,R), which can be looked up using the command clidata in CLIFFORD. 
The main reason for using these isomorphisms is that R ® Ci P:q ~ Ct p q , so we 
don't have to worry about K tensor products of spinor modules and use only 

S® M S* ~Mat(2 A ',R). 

Example 2. A spinor basis for Ct\ \ ~ Mat(2,R) used in the isomorphism 6), given 
the orthogonal generators e\ = 1 = — e% for Ci\ : \ and a primitive idempotent /, 
maybe chosen as 

S = Ct XA =span R {/ = i(l +ei, 2 ),ci/= +e 2 )}- (10) 
In this basis, the following matrices represent the basis elements in Ct\ \ : 




In general, the images of the additional generators e, ®e\e2 needed to generate 
Cl P: q(g)C£i,i ~Mat(2,Cf Pii:/ ) are given by ei®\e\e2\ with entries from Cl PA , and read 

e ( - ® eie2 — ^ . The isomorphism 6) gives, via iteration, Ct p . q ~ Cl p - q $ ® 

Mat(2«,R) if p > q and Cl PA ~ CX , <7 _ p (8)Mat(2P,R) if 9 > p. 



3 Computing with CLIFFORD and Bigebrain tensor algebras 

We provide examples of Maple code how to set up tensor products of Clifford alge- 
bras in CLIFFORD and Bigebra. For the usage of these packages see UBS), the 
help pages which come with the package, and the package website [3 |. The Maple 
worksheets with code for the described methods will be posted at |6|. 

Loading the package using >with (Clifford) ; with (Bigebra) ; exposes the 
exported functions. To set up a Clifford algebra, say Cli.i, one needs to define the di- 
mension >dim_v:=2+2; and the bilinear form >b : =iinaig [diag] (l$2,-l$2);0 
Basis elements ej are written as strings, e.g., eiwe4 stands for e\ Ae4, etc., whereas 
id stands for the identity of the Clifford algebra. Bigebra exports also the 
(graded) tensor product &t, which is multilinear and associative. Then, a tensor 
product ei 2 ®<?i reads &t (eiwe2 , el ) , and permutations of tensors are implemented 
by maps >switch(&t (ei,e2) , l) = &t(e2,ei) (the ungraded switch) or, in the 
graded case, >gswitch (st (el, e2) , l) = -st(e2,el) (the graded switch). The 
extra index ; in either switch (here we have used l in each) tells [g] switch to 
swap the i-th and the (i+l)st elements. Again, you get help by typing >?switch and 
>?gswitch at the Maple prompt. 



The default name of the bilinear form in CLIFFORD and Bigebra is B, however other names 
can also be used. So, when the bilinear form B is left undefined (unassigned), computations are 
performed in a Clifford algebra d(S) for an arbitrary bilinear form B (see, e.g., (5||13I ). 
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The Clifford product cmul by default implicitly uses the bilinear form B as in, 
for example, >cmui (el, e2 ) =eiwe2+B [ i, 2 ] *id. However, it can also use B or 
any other Maple name explicitly as an optional argument, e.g., >cmui [K] (el, e2) 
= eiwe2+K[i,2]*id, allowing to compute in different Clifford algebras in the 
same worksheet. Let B, B\, B 2 hold the bilinear formal for C£ p+l ^ +s ,C£ Pt q and 
Clr. s , and let bas2GTbas be the graded algebra isomorphism (0) given by ej G 
d P ,q i-^&t (ei, id) (/ C {l,...,p + q}) and ej G C£ r ,s (id,eJ) (J C {l,...r + 
s}), then the following procedure 



cmulGTensor : 


= proc (x, y, Bl, B2) 




local 14; 






f4:=(a,b, 


x, y) ->cmul [Bl ] (a,b) , cmul [B2; 


(x,y) : 


eval (subs 


( '&t '=f4, gs witch (St (x, y) ,2) 


) ) ; 


end proc : 







implements the graded tensor product algebra as explained in section |2~41 We prove 
in the worksheet cmulGTensor .mw Q that 



bas2TGbas (cmulGTensor (X, Y, Bl, B2) ) = 

cmul [B] (bas2GTBas (X) ,bas2GTbas (Y) ) 

is the graded algebra isomorphism (|4) with the inverse GTbas2bas. The grade and 
reversion involutions work as expected. The limit 6-)0 implements the wedge 
product on A V = f\V\ ® /\V 2 . 

We discuss in more detail the ungraded tensor product for C£ p+ \^ + \ ~ C£ p , q (S 
Cl\\, the isomorphism 6) of Thm. [2] which we call bas2Tbas while its inverse is 
Tbas2bas. We have in Ct\ \, with generators e\ = 1 = — e|, the volume element 
ft> = e i e 2 with ft) 2 = A = 1, as we have to use the bilinear form jQi of Thm. [3] 
so we can still use jB\ = B\. The isomorphism 6) reads e/ G C£ p . q >— Ht (ei,w) 
(7C {l,...,p + q}) andey G C£ M i-^&t (id,eJ) (J C {l,2})@The tensor Clifford 
product is graded isomorphic to the Clifford product on C(, p+ \ xl+ \. For a proof see 
the worksheet cmuiTensorii .mw. (6 1 The following procedure implements the un- 
graded tensor product algebra as explained in section l2~5l 

cmulTensor : = proc (x,y,lBl,B2) 
local f4; 

f4: = (a,b,x,y) ->cmul [1B1 ] (a,b) , cmul [B2] (x, y) : 
eval ( subs (' St ,= f 4 , switch (&t (x, y) , 2 ))) ; # ordinary switch 
end proc : 

The isomorphism bas2Tbas is this time more involved, as are the grade and rever- 
sion involutions. We still get the isomorphism 

bas2Tbas (cmulTensor (X,Y,B1,B2) )= 

cmul [B] (bas2TBas (X) , bas2Tbas (Y) ) 

proved by direct computation explicitly. Further details are provided in the work- 
sheet. Note that we do not have the limit B — > 0, as Qi ~ B 2 needs to be nonde- 
generate, and the naive limit replacing cmul by wedge gives false results. In the 

5 Computations in the worksheet cmulGTensor .mw are performed for arbitrary not necessarily 
symmetric or diagonal bilinear forms. 

6 Here, w in St (el, w) stands for the volume element co. Then, in the code of cmulTensor, the 
form \B\ is denoted as 1B1. 
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worksheet we show how to produce the Grassmann basis for the tensor algebra, and 
how the isomorphism and the involutions work. 

In each worksheet we have benchmarked computations using the generic cmul 
routine for the Clifford product from CLIFFORD (in dim V < 9) versus the Bigebra 
tensor routines cmul [G] Tensor. For orthonormal bases we roughly get equal run- 
times. The more complicated data structures of the tensor algebras is compensated 
by not computing some of the off-diagonal terms. 

We add to this brief discussion that it is easily possible to iterate this morphism, 
and to provide CLIFFORD code for computations in Clifford algebras 

a p+Kq+k ~ ci M ® a u ® • • • ® a hh 

" V ' 

k factors 

or use the mod 8 periodicity. 



4 Computations using matrix algebras over Clifford numbers 

The isomorphism 6) from Thm. 2 was explicitly defined by Lounesto in |fl6l 
Sect. 16.3]. We will use this matrix approach to perform computations in Ct%2 — 
Mat(2,Q?7i) Q. Let {e\,...,e%} be an orthonormal basis of R ' generating the 
Clifford algebra Cijj such that ef = 1 for 1 < i < 7, e| = —1, and e,e ; = — eye,- for 
< i,j < 8 and ; ^ j. The following 2x2 matrices (compare with (fill ) 

£ '=(o'4) for ' 1 8 - £9 =(ii)< Eio =(?"i) (i2) 

anti-commute and generate C£g.2 Qln order to effectively compute in (3?8,2, we define 
the isomorphism <p : Mat(2,(3?7.i) — > Cl%2 as a Maple procedure phi from the asvd 
package0 Its inverse is given by the Maple procedure evaim. In the first step, we 
compute 1,024 2 x 2 matrices £> with entries in Cl-i.\ which represent the basis 
monomials e\ = e\ x ei 2 ■ ■ -ei k , i\ < ii < ■ ■ ■ < ik, < k < 8. We store these matrices 
in a list 38. 

For example, for the basis 3s ^ of A:-vectors in CLIFFORD we compute all [ k J 
products Ei { &cm£/ 2 scm • • • son Ej k where &cm is a CLIFFORD procedure to com- 
pute a product of Clifford algebra-valued matrices with the Clifford product applied 
to the matrix entries. Once the matrix representation (Q~2} has been chosen@the list 
3§ ~\J k 3§k can be saved and read into a next Maple session thus avoiding the need 
to repeat this step. 



7 Notice that Ef = Eg = 1 for 1 < i < 7 and £| = £^ = — 1 in this explicit representation. Also, 
0?i : i ~ Mat(2,]R) is such that the matrices e\ = ^ J J and ei = ^ generate C£i : i. 

8 The asvd package is part of the CLIFFORD library. It was introduced in (2). 

9 Note that other representations are possible. See II 1 61 Sect. 16.3]. 
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Having computed the basis matrices 38, we are now ready to compute in (3?8,2- 
For example^ let x = 2ID + 4E U23 - 10£i,5, 7 ,8,io and y = ID + 4£ 1 , 2 ,3.7 + 
£1,5,6.8 — 3£i 467. We can find the Clifford productx &CMj of x and y in Cl%^ using 
the following procedure: 

' &CM' : =proc (x : : algebraic , y : : algebraic ) 
local xy; global phi , BBB ; 

if not type (evalm (x) , climatrix) then error 

'evalm (x) must be of type climatrix' end if: 

if not type (evalm (y) , climatrix) then error 

'evalm (y) must be of type climatrix' end if: 

xy : = displayid (evalm (x) scm evalm (y) ) ; 

return phi (xy, BBB); 
end proc : 

Once the basis S3 is stored in Maple as a list bbb, the procedure &cm treats it as a 
global variable and it return sT*! 

— 40£ , 23 : 5 : 8 : 10 + 2£'i,5,6,8 - I6.E7 — 4Ei j2 ,3 + 10^1,5,7,8,10 — 10^6,7,10 

+ 8£l. 2.3,7 + 4&2,3,5,6,8 — 12^2,3,4,6,7 — 30£4.5.6.8,10 — 6£[ 4.6.7 — 2ID . (13) 

Details of the above computations can be found in the worksheet G82 .mws J6|. In 
the worksheet we also define the respective grade and reversion involutions, and the 
graded algebra isomorphisms Cl%2 Mat(2,Q?7 i) are given by phi and evalm. 
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